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Abstract 

We address the question whether Bohmian trajectories exist for ah times. 
Bohmian trajectories are solutions of an ordinary differential equation involv- 
ing a wavefunction obeying either the Schrodinger or the Dirac equation. Some 
f^ ' trajectories may end in finite time, for example by running into a node of the 

^J^ . wavefunction, where the law of motion is ill-defined. The aim is to show, under 

vQ I suitable assumptions on the initial wavefunction and the potential, global exis- 

<^ ' fence of almost all solutions. We provide a simpler and more transparent proof 

^^ . of the known global existence result for spinless Schrodinger particles and extend 

the result to particles with spin, to the presence of magnetic fields, and to Dirac 
wavef unctions. Our main new result are conditions on the current vector field on 
configuration-space-time which are sufficient for almost-sure global existence. 
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1 Introduction 

We study a mathematical question arising from and relevant to Bohmian mechanics 
[SI HI HH ini HHj and its variant based on the Dirac equation jHl IZ| (henceforth referred 
to as the "Bohm-Dirac theory"). In these theories, the motion of particles is defined 
by ordinary differential equations (ODEs) involving the wavefunction, see Q and (0) 
below. The mathematical question we address is global existence, i.e., whether (under 
what conditions and how often) the particle trajectories are well defined for all times. 
One obstruction to global existence is that the velocity given by (jH)) or (j3)) is singular at 
the nodes (i.e., zeros) of the wavefunction. In particular, there are trajectories that are 
not defined for all times because they run into a node. Thus, the strongest statement 
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one can expect to be true is that global existence holds for almost all solutions of the 
equation of motion. As we show, this is in fact true for suitable potentials and initial 
wavef unctions. As a by-product, one obtains from almost-sure global existence the 
equivariance of the |^ip distributions. 

The relevance of Bohmian mechanics to the foundations of quantum mechanics arises 
from the fact that a world governed by Bohmian mechanics satisfies all probability rules 
of quantum mechanics HI El IHl OH] • Bohmian mechanics thus provides an example of 
a "quantum theory without observers," one in which no reference to observers is needed 
for the formulation of the theory, and an explanation of the quantum probabilities in 
terms of objective events. 

The authoritative paper on global existence of Bohmian trajectories is by Berndl et 
al. |lj; see also |2]. We note that the proof given by Holland ^3 p. 85] is incorrect 
(see |1] for details). We also remark that the general existence theory for first order 
ODEs with velocity vector fields that are not Lipschitz but only in some Sobolev space 
P] does not apply to Bohmian trajectories. The results of [H] hold for vector fields 
with bounded divergence, while the divergence of a Bohmian velocity field, such as in 
Q and (JH), typically diverges at nodes of the wave function. Berndl et al. jl] already 
proved almost-sure global existence for suitable potentials and initial wavef unctions; 
while they give a proof only for spinless nonrelativistic particles, a similar proof could 
presumably be devised for Bohmian mechanics with spin ^1 Ej and the Bohm-Dirac 
theory. We provide here an alternative proof that is shorter and more transparent than 
the proof by Berndl et al. Our result covers all cases covered by their existence theorem; 
in addition, our result also covers Bohmian mechanics with spin and magnetic fields 
and Bohm-Dirac theory; for the latter our result and its proof become particularly 
simple thanks to the fact that the Bohm-Dirac velocities are bounded by the speed of 
light. Even more generally, our result can be applied to any Bohm-type dynamics, as 
we formulate conditions on the current vector field on configuration-space-time that are 
sufficient for almost-sure global existence. 

There are three ways in which a trajectory can fail to exist globally: it can approach a 
node of the wavefunction (where the equation of motion is not defined) , it can approach 
a singularity of the potential (where the equation of motion need not be defined), or it 
can escape to infinity in finite time. Hence, the main work of any existence proof for 
Bohm-type dynamics is to show that almost every trajectory avoids the "bad points" 
(nodes, singularities, infinity) in configuration space. The method of Berndl et al. is 
based on estimating the probability fiux across surfaces surrounding the bad points and 
pushing these surfaces closer to the bad points; in the limit in which the surfaces reach 
the bad points, the fiux vanishes. 

The advantage of our approach is that it does not require skillful estimates and does 
not involve limits. Instead, our method is based on considering a suitable nonnegative 
quantity along the trajectory that becomes infinite when the trajectory approaches a 
bad point; if such a quantity has finite expectation, at least locally, then the set of initial 
configurations for which it becomes infinite must be a null set. That the expectation be 
locally finite can be paraphrased as an integral condition on the current vector field. 

To illustrate our method, we briefiy describe an argument of this kind: the total 
distance D traveled in configuration space in the time interval [0,T] becomes infinite 



when the trajectory escapes to infinity during [0,T]. To prove that D is almost surely 
finite, we prove that it has finite expectation. A calculation shows that 
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where J is the spatial component of the current vector field. Thus, the finiteness of 
the right hand side of (0) is a natural condition on the current ensuring that almost no 
trajectory escapes to infinity in [0,T]. 

This argument was already sketched in ^3] ; it was inspired by a similar consideration 
in the global existence proof of p[ll for Bell's jump process for lattice quantum field 
theory, another Markov process depending on a wavefunction ipt and having distribution 
|-^tp at any time t; the quantity considered there was the number of jumps during [0, T]. 
Finally, a related argument was also described in Remark 3.4.6 of ^, see our Remark 5 
for a comparison. 

Our method is also in a way more elementary than that of Berndl et al.: we do 
not make use of the nontrivial fact that the iV'p-probability of crossing a surface S in 
configuration-space-time is bounded by Jj^lda ■ j\ where da is the normal on S with 
length equal to the area of the surface element and j is the current vector field. Indeed, 
we use this fact only for surfaces lying in t = const, slices of configuration-space-time, 
for which it is much simpler to prove, see Lemmas To be sure, the statement about 
general surfaces is interesting in its own right and also relevant to other applications 
such as scattering theory, but its proof takes several pages in j^j. 

While our innovation concerns sufficient conditions on the current for almost-sure 
global existence, there remains the functional analytic question of deriving these condi- 
tions from conditions on the potential and the initial wavefunction. We carry this out in 
several example cases but contribute nothing original; we employ the same arguments 
as Berndl et al. or standard arguments. 

This article is organized as follows. In Section |21 we give the definition of Bohmian 
trajectories for both the Schrodinger and the Dirac equation; we elucidate the relevance 
of the current vector field to the trajectories and their distribution. In Section El we 
state and prove our results in terms of a current vector field. In Section 0] we state and 
prove our results for Bohm-Dirac theory. Finally, in Sectional we state and prove our 
results for Bohmian mechanics. 



2 Setup 

We briefiy recall Bohmian mechanics and the Bohm-Dirac theory for a system of A^ 
particles. Then we describe what singularities we will allow in the potential. Finally, 
we point out how for both Bohmian mechanics and Bohm-Dirac theory the trajectories 
arise from a current vector field on configuration-space-time. 



2.1 Equations of Motion 

In Bohmian mechanics, the wavefunction is a function ip : M. x M.^^ -^ C^ where M 
represents the time axis, M^^ the configuration space of N particles, and C'^ the value 
space of the wavefunction representing the internal degrees of freedom of the particles 
such as spin (and possibly quark fiavor etc.). ip = ip{t, q^, . . . , q^) evolves according to 
the Schrodinger equation 

where mi and Cj denote mass and charge of the z-th particle, c the speed of light, A 
is the external electromagnetic vector potential, and V is the potential, which may be 
Hermitian k x fc-matrix valued. For particles with spin in the presence of magnetic 
fields, the potential includes a term ^^ 2m^(^ ^ ^)(9j) ' '^i where ctj is the vector of 
spin operators (Pauli matrices for spin-|) acting on the spin index of particle i; this 
form of the Schrodinger equation is known as the Pauli equation. 

The law of motion for the trajectory Qi{t) of the i-th particle reads 

at rrii ip* ip 

where Q = {Qi, ■ ■ ■ , Q^) is the configuration, and ip*(t) denotes the inner product in C^. 
The right hand side of is ill-defined when and only when either 'il){t, Q) = Q (node of 
ip) or ip is not differentiable at (t, Q). For an explicit example of a trajectory that runs 
into a node of ip, see jl]. 

In Bohm-Dirac theory, the wavefunction is a function ^ : M x R'^^ ^ C^ = (C^)'^'^ 
evolving according to the Dirac equation 

ih— = - ^ ichcxi ■ Vq^ip + 1/(^1, . . . , q^)ip, (4) 

1=1 

where cxi denotes the vector of Dirac alpha matrices acting on the spin index of particle z; 
we have included the mass terms in the potential V , which is Hermitian 4^ x 4^-matrix 
valued. In the presence of magnetic fields, V includes a term — ^^ eiA{q^) ■ oti. 
The law of motion for the trajectory Qj(t) of the z-th particle reads 

The right hand side is ill-defined at nodes of ip and only there. 

2.2 Singularities of the Potential 

Among the physically relevant examples of potentials V = V{qi, . . . , g^) is the Coulomb 
potential, 

Viq„...,q^) = J2j-^ (6) 



which is singular at coincidence configurations (those with g^ = q, for some i ^ j). 
This motivates us to allow that V is defined only on a subset Q C M^^; e.g. in the 
case of Coulomb interaction, Q is the set of non-coincidence configurations; in the 
case of an external Coulomb potential generated by charges located at 2:1, ... , zm, Q = 
(M^\{zi, . . . , zm})^ ■ One cannot expect a Schrodinger wavefunction to be differentiable 
on the singular set M'^^ \ Q of the potential, as exemplified by the ground state of the 
hydrogen atom, which is proportional to exp(— A|qf|) for suitable A > 0. Thus, the right 
hand side of (jH)) may be ill-defined on M^^ \ Q, and we will use differentiability of ip only 
on Q. For the Coulomb interaction and the external Coulomb potential, Q is of the 
form Q = M"^ \ Uj^iSg, where Si are hyperplanes. Our method of proof allows somewhat 
weaker assumptions: 

A closed set S* C M*^ is admissible, if there is a 5 > such that the distance function 
q I— »■ dist(g, S) is differentiable on the open set (S* + 5) \ S, where S + 6 = {q E M.'^ : 
dist(g, 5*) < S}. Then the configuration space Q is 

m 

either Q = M%r Q = M"^ \ |J ^^, (7) 

e=i 

where Si, . . . , Sm are admissible sets. For example hyperplanes are obviously admissible 

sets. 

2.3 The Current Vector Field 

There is a common structure behind the laws of motion (jH} and (jSj): they are of the 
form 

dQ... _ J{t,Q{t)) 

dt^^ f{t,Q{t)) ^'^ 

where j = (j°, J) is the current vector field on configuration-space-time M x Q, defined 

by 

J = (iV^p, ^Jmr{V,^ - t A(q,))V^, . . . , ^lmr{V,, - ^A(q^))V^) (9) 
in the Schrodinger case and 
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,cilj*aitlj,. . . ,cilj*a.Ntlj] (10) 

in the Dirac case. Provided that ip is sufficiently differentiable, j has the following 
properties, which we take to be the defining properties of a current vector field: 

j = (jO, J) is a C^ vector field on M x Q (11a) 

d 
divj = 5^9^j^ = (lib) 

At=0 

j° > whenever j 7^ (He) 

dqf{t,q) = l \/teR. (lid) 

Q 



We will call points in ,jY = {(t,q) ^ M. x Q : j(t,q) = 0} the nodes of j. We write 
^ = {q & Q : j(t, q) = 0} for the set of nodes at time t. 

Let Qq{t) denote the maximal solution of (jH)) starting in g G Q\ jVq defined for 
t G {T^iTq). It is a reparameterization of an integral curve of j, see Remark |H1 for 
details. We will formulate our existence theorem first purely in terms of the current 
vector field, and then apply our result to the currents and (fTUI) . 

2.4 Equivariance 

We now explain the notion of equivariance, and what needs to be shown to prove equiv- 
ariance. We first remark that equivariance is a crucial property of Bohm-type dynamics, 
in fact the basis of the statistical analysis of Bohmian mechanics [llj and thus the basis 
of the agreement between the predictions of Bohmian mechanics and the prescriptions 
of quantum mechanics. We also remark that, while full equivariance will be a conse- 
quence of the existence result, a kind of partial equivariance can be obtained before, see 
Lemma n below; our existence proof will exploit this partial equivariance. 

Before we define equivariance, we introduce some notation. Let £^{Q) denote the 
Borel cr-algebra of Q. Let /ij be the measure on s^{Q) with density j°(t) relative to 
Lebesgue measure, 

fj.t{B) = jdqfit,q) (12) 

B 



for all B G ^(Q). By (jlld|l . yUj is a probability measure; in Bohmian mechanics and 
Bohm-Dirac theory, fit is the |'?/'(t)P distribution. We introduce a formal cemetery 
configuration ^ and set Qq(t) := <:^ for all t ^ {r^ , r+), respectively, if (0, q) is a node of 
J) Qq(t) '■= C> for all t 7^ 0. Let v^t : Q — > Q U {<C>}, ^til) = Qq{t), denote the flow map 
of (jHl), and letyjiMx Q^Mx(QU {(}}) be the flow map on configuration-space-time 
defined by ip{t,q) = (t,<^i(g)). Let Qt = {q E Q\ ^o ■ r~ < t < r+} = ip^\Q). 
Standard theorems (see, e.g.. Chapter II of [inj) on ODEs imply that ip is C^ on the 
maximal domain {(t,q) G M x (Q \ ^) : ipt{q) 7^ (^}, which is open; in particular, also 
Qt is an open set. 

Let pt be the distribution of Qq{t) if q has distribution /xq, i-e., 

pt = /io o ip~^ . (13) 

One says that the family of measures fit is equivariant on the time interval I ii pt = Pt 
for all t E I. (The interval I may be finite or infinite.) 

Let J^t '■= ^t{,Qt) = ft{Q) n Q be the image of the flow map in Q at time t. The 
following lemma formulates "partial equivariance." 

Lemma 1 Let j = {f,J) satisfy ffTTal . (jllbj) . and fTTc|l . Then for all B G £/{Q) and 
all t G R, 

pt{B) = pt{B n Jt) ■ (14) 

We know of two ways of proving this lemma, requiring comparable effort. One 
proof, given in ^ and in more detail in [2j, goes as follows, pt has a density that obeys 
a continuity equation, and j° satisfies the same continuity equation. By uniqueness of 
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solutions of this linear partial differential equation, one obtains that j'^it) coincides with 
the density of pt on J^f . An alternative proof, which we give below, is based on applying 
the Ostrogradski-Gauss integral formula to j on a cylinder formed by the trajectories 
over a polyhedron in Q. 

Proof of Lemma Ql Without loss of generality, t > 0. For any (i-chain of singular 
simplices E in Qt, the cylinder F formed by the trajectories over E, F = (f{[0,t] x E), 
is a d + 1-chain in configuration-space-time M x Q. Applying the Ostrogradski-Gauss 
integral formula to j and F, we obtain 



= / dtdqdivj 



da- j 



OF 



where da is the outward pointing surface normal with length \da\ equal to the area of 
the surface element. The surface dF of the cylinder consists of three parts: the mantle 
(/9([0,t] X dE), the hd f{{t} x E), and the bottom (f{{0} x E). The integral over the 
mantle vanishes as the mantle consists of integral curves of j and is thus tangent to j. 
The integral over the lid is / ,^. dqj^{t, q) and that over the bottom is — J^dqj^{0,q). 
Therefore, we obtain 

= ptME)) - po{E) = pM{E)) - ptME)) . 

Any two measures that agree on the d-chains (and thus in particular on the compact 
rectangles) agree on a fl-stable generator of the cr-algebra s^{Qt) and are, by a standard 
theorem, equal. Since v?t is a bijection Qt -^ J't-, we obtain (fT^ . D 

What remains to be shown to prove equivariance is that pt{Q\ •^t) = 0. 



3 A General Existence Theorem 

Let e^(Q) denote the set of all bounded Borel sets in Q. 

Theorem 1 Let Q C M.'^ be a configuration space as defined in (|7j) and let j = {j^, J) 
be a current as defined in ()11|). Let T > and let (ft '■ Q -^ Q^ {(}} denote the flow 
map of (jHl). Suppose that 



d_ J_ 



'^,)fit,q) 



\/B e ^(Q) : dt dq 

MB)\{0} 

T 

WBE^iQ): f dt f dq 
MBMO} 
and, if Q = W^ \ UiSg, in addition that for every £ G {!,..., m}, 



< CXD 



Jit, q) ■ n 

m 



< cxo - 



(15) 



(16) 



3(5 > V5 G =^(Q) : I dt I dql(qe{Se + S))^-^^%4^^^<oo. (17) 

' ' dist(g, S^^) 



MB)\{0} 



Here dist(g, S() is the Euclidean distance of q from Si and ei{q) = — Vgdist(g, 5"^) is the 
radial unit vector towards S^ at q & Q. Recall that for 6 sufficiently small the distance 
function is differentiahle on Si + 6. 

Then for almost every q E Q relative to the measure fio{dq) = j^{0,q)dq, the solu- 
tion of (jHl) starting at Q{0) = q exists at least up to time T, and the family of measures 
f't{dQ) = j^{'tj(l)dq is equivariant on [0,T]. In particular, if fllSj) . fll6|) and, if appropri- 
ate, (fTTj) are true for every T > 0, then for fiQ-almost every q E Q the solution of (jH)) 
starting at q exists for all times t > 0. 

Remarks: 

1. We can formulate the meaning of each of the conditions (fT3j) . fll6p . and (fT7|) as 
follows. If (|T5|l holds, then /io-almost no trajectory approaches a node during 
[0,T]. If (lTI)jl holds, then /iQ-almost no trajectory escapes to oo during [0,T]. If 
(fTTjl holds, then /iQ-almost no trajectory approaches a point in the singular set 
Wp^-^Se during [0,T]. 

2. To obtain existence also for negative times, one can apply Theorem d to the time 
reversed current 

jit,q)={fi-t,q),~Ji-t,q)). (18) 

The integral curves of J are the time reverses of the integral curves of j. Obviously, 
with j also J satisfies (fTTj) . If J satisfies (fT3jl . (fT^ . and, if appropriate, (fTTj) for 
T > 0, we obtain almost-sure global existence of Qg(t) on [— T, 0]. 

3. It suffices to consider in (fT3j) . (fTB|) and (fT7|) for the sets B instead of all bounded 
Borel sets just the balls around the origin. This is because enlarging B cannot 
shrink the integral. For the same reason, it suffices to integrate over Q\^ instead 
of the not easily accessible sets ft{B) \ {(}}■ 

4. Actually the proof of Theorem Q] works in the same way with the following slightly 
weaker conditions. Instead of (fT^ it suffices to assume that 



WB e^{Q) 3R<oo : f dt f dq l{\q\ > R) 

MB)\{0} 

and p7|l can be replaced by 






< oo 



WBe^iQ) 3S>0: [dt [ dql(qe{Se + 6))^-^^^j^4^^0^<oo. 

J J dist{q, Si) 

MB)\{0} 

We chose to state the theorem with the stronger assumption to simplify the pre- 
sentation and because the weaker assumptions will not be used in the following. 

Proof of TheoremUi Let Qq{t) be the maximal solution of (jH)) starting in q, as described 

in Section 12.31 Since we deal only with positive times in the following, we write Tq for 

r+. 
q 



There are three ways in which Qq(t) can fail to exist globally: the trajectory can 
approach a node, approach a point on the singular set USg, or escape to infinity in 
finite time. More precisely, if q ^ ^ and Tg < oo there exists an increasing sequence 
(^n)ngN with tn — * Tg such that either there is x & jVr^ U \Si=\^i. with Qq{tn) -^ x or 

\Qq{tn)\ -^ OO. 

To see this, suppose that Tg < oo and that such a sequence did not exist. Then 
Qq '■= {(t,Qg(t)) : t G {0,Tq)} C (M X Q) \ ^ would remain bounded and bounded 
away from the complement (M x USe) U ^/K. Since (M x Q) \ ^ is open, there would 
be a compact set ii' C (M x Q) \ ^ such that Qg C K°, with K° the interior of K. 
However, the vector field (1, J/j^) is C^ on (M x Q)\ ^ and thus uniformly Lipschitz 
on K. Therefore, all of its maximal integral curves either exist for all times or hit the 
boundary of K, in contradiction to the hypotheses. 

Let now q ^ ^ and Tq < T. If there is x G ^^ and (t„) such that Qq(tn) — ^ x, 
then j^(tn, Qqitn)) -^ 0. Hence, the total variation of t i— > log j°(t, Qq(t)) up to time T 
diverges, i.e., Lq = oo where 

min(T,T,) 

Lq= j rft|^logj°(t,Q,(t))| forgGQ\^o. (19) 



We now show that L^o '■= {(1 & Q\ -^q : L^ = oo} is a /iQ-null set. For this it suffices to 
show that for any bounded set B G e^(Q), B fl L^o is a /io-nuU set. For this in turn, it 
suffices that the average of Lq over B relative to the measure /^o be finite: 

mm(T,Tq) 
B BO 

[the order of integration can be changed since the integrand is nonnegative] 



T 



dt / dqf{0,q)l{Tg > t) 



d aO 



3\t.QM\ 



f{t,Qg{t)) 

Ud/dt + j,-Vq,)f{t,q')\ 



= Jdt J dq' \id/dt + j, ■ Vq>)f{t,q')\ 9 oo. 

MB)\{0} 

This shows /io(-Z^oo) = and thus that the solution Qq(t) of (jHI) /io-almost surely does 
not approach a node during [0,T]. 

Now we consider the cases that either 





B 




T 




= dt / pt{dq 




MB)\{0} 


by Lemma H] 





lim \Qq(tr, 



OO . 



or 



3xGU^iS^: QgiQ^x. 



Hence, for such initial conditions either the total variation of t i-^ IQg(^)l is infinite, i.e., 

mm{T,Tg) 



Dg = oo where 



D, 



dt 



j,\QM 



for g e Q \ ^ , 



or the total variation oi t \-^ logdist((5(j(t), Se) restricted to Se + 6 is infinite for some 
£ G {1, . . . , m} and any 5 > 0, in particular for the one in (|T7jl . i.e., Vq^i = cxd where 



min(T,r,) 



v„. 



I dtl{Qg{t)e{Se + 6))\j^\ogdist{Qg{t),Se)\ for g G Q \ ^ . 



Therefore it suffices to show that D^o :={(?€ Q \ ^ : Dg = oo} and \4o/ := {q ^ 
Q \ ^ : Vg^i = 00} are /iQ-nuU sets, and for this we proceed as for L^o- 

Let B G e^(Q). Then (jHJ, followed by exactly the same arguments as for Lg, shows 
that local expectations of Dg are finite, i.e. 



dqfiO,q)Dg 



q' 






< CXD 



dt / dq' 

B MB)\{0} 

Hence /io(-Doo) = 0. For local expectations of Vg ^ we obtain, again with (jS)) and LemmaC] 



dq3'{Q,q)Vg, 



min(r,Tg) 

dqf{Q,q) [ dtl{Qg{t)eiSe + 6)) 



B 



B 







Qg{t) ■ eeiQgit)) 



dist{Qg{t),St> 



dt / dq'l{q' e {Se + 6)) 
MB)\{0} 



J{t,q')-e,{q') 



dist(g'. Si 



ITtI 
< 00 . 



Hence also /Uo(K»/) = 0, concluding the existence part of the theorem. 

It remains to show equivariance. Since the probability of reaching <) before time T 
vanishes, we have pt{Q) = 1 for all t G [0, T]. Since pt < pt by Lemma [Hand pt{Q) = 1 
by (jlldj) . we must have pt = pt, which is equivariance. D 

Remarks. 

5. A reasoning closely related to our method of proof is also applied in ^, Re- 
mark 3.4.6. There, an expression analogous to p5|l is used to control the proba- 
bility of reaching an e-neighborhood of ^ before letting e -^ 0. Apart from the 
fact that the argument is applied there only to the nodes and not to singularities 
and infinity, it is also unnecessarily complicated, mainly because it considers an 
e-neighborhood instead of fully exploiting the integral p9|) . 
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6. The proof of equivariance was the only place where we used the property (jllc 
of a current vector field. The existence statement of Theorem ^ holds as well if j 
satisfies (fTT|) except for ()lld|) : in particular, we may allow fit{Q) = oo. 

7. Here is another equivariance result that does not use (jlldjl : Let Q be a con- 
figuration space as in ((Zj) and let j satisfy f|TT|l except for fllldjl . Suppose that 
almost-sure global existence holds in both time directions, starting from any time. 
Then the family of measures fit is equivariant on M. 

To see this, note that for equivariance we need to show merely that Q \ J^^ is a 
/i(-null set, or, in other words, that for /ij-almost every g G Q the integral curve 
of j starting in (t, q) reaches back in time to time 0. But this is immediate from 
almost-sure global existence in the other time direction, starting at time t. 

Thus, if both j and J as defined in (J18p and their time translates satisfy (jl5|) . (J16p . 
and, if appropriate, (jl7|) for all T > 0, we obtain equivariance without ()lld|) . 



8. Condition (fT6|) can be replaced by the condition 

the first order derivatives of J are bounded on [0, T] x Q. (20) 

To show this, we show that under this assumption every unbounded solution Qq{t) 
with Tq <T has Lg = oo, with Lg defined in (ITTHl . 

To see this, first note that the solutions of (jH)) are reparameterizations of the 
integral curves of j. In more detail, let 7g(s) = (7°(s), rq(s)) be the unique 
maximal integral curve to j, 

'-^=m^)). (21) 

starting in (0, g) G M x (Q\^) and defined for s G {cr~ , a^). Since j° > outside 
nodes, 7g(s) is monotonically increasing, and hence the map 

s s 

r d'y^ r 



t,is) = 7,"(.) = I d~s^= I dsfi^gis)) 

is invertible on its image {T~,r^), where r^ = lim^^^± tg(s), with inverse Sq{t). 
Since 

Qq(t) = Tq{sq{t)) IS the uulque maximal solution of (jHl) with Qq{0) = q; it is 
defined for t G (r~, r+). 

Now suppose that \Qq(tn)\ — > cxd for some tn — >■ t+. Then also \Tq{sq{tn))\ — >■ cxo. 
Since the derivatives of J are bounded, there are constants A, R > such that 
\J{t,q)\ < A\q\ for all t G [0,T] and all g G Q with |g| > R. Since dTq/ds = 
J(7q(s)), it follows that |rq(s)| < max(|rg(0)|, i?) e^''; thus, an integral curve of j 
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cannot escape to spatial infinity in a finite interval of the parameter s; in other 



words, cr+ = oo. But then 

POO 

< = / dsf{-fqis)) < CX) 
Jo 

implies the existence of an increasing sequence (s„) with s„ -^ cxo such that 
j°(7g(s„)) — s> 0, and therefore Lg = oo. 

4 Global Existence of Bohm— Dirac Theory 

The Dirac Hamiltonian for N particles is 

N 

where we assume a nonsingular V G C°°(M^^, Herm(C^ )). According to ^, Hd is 
essentially self-adjoint on C^(]R^^,C'' ) and we denote by Hd the unique self-adjoint 
extension. 

Since the Dirac matrices a have eigenvalues ±1, the velocities in (0) are bounded 
by c. Consequently, the Dirac current ()10j) satisfies \J\ < cyN j^. This fact makes the 
proof of global existence particularly simple, as expressed in the following corollary to 
Theorem H 

Corollary 1 Let Q = W^ and let j = (j°, J) be a current as defined in (jllj) . Suppose 
that there is a global bound on velocities, i.e., a constant c > such that \J\ < cj^. 

Then for fiQ-almost all q G W^ , the solution of (jHl) starting at (5(0) = q exists for all 
times, and the family of measures fit is equivariant. 

Proof of Corollary We show that assumptions p5|) and p6p of Theorem Q] are 
satisfied for any T > 0. The key observation is that due to the bound on velocities, 
bounded sets in configuration space stay bounded under the flow. More explicitly, for 
any bounded set B G J^{M.'^) contained in, say, the ball B^ of radius r > around the 
origin, (pt{B) \ {(}} will be contained in Br-^.ct and thus in Br+cT provided t G [0,T]. 
Now \{dt + 4j ■ Vg)j°| < |9ij°| + c|Vgj°|, and the functions |J|, |9ij°|, and c|Vgj°| 
are continuous and therefore bounded on the compact set [0, T] x Br+cT- Hence the 
integrals in (fT3j) and ()16|) are finite. This implies existence for all positive times. For 
negative times apply the same argument to the time-reversed current J, for which the 
same velocity bound holds. D 

Applying Corollary Q to Bohm-Dirac theory, we obtain global existence of Bohm- 
Dirac trajectories under very general conditions. 

Theorem 2 Let V e C°°{M.^^, Herm(C^'^)) and ip{t) = e-'^^^ipiO) with ^(0) G 
C°°(M3^,C4'^)nL2(R3A^,C4'') and ||^(0)|| = 1. 

Then the solution Qq(t) = {Qi(t), . . . , QAr(t)) of © with Qq{0) = q exists globally 
in time for almost all q G M^^ relative to the measure fio{dq) = \4'{0,q)\'^dq, and the 
\4'{t)\'^ distributions are equivariant. 
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Proof. According to [HI, for ^(0) e C^ (M.^^ , C^"" ) one has z/'(t) G C^(R^^,C^^) and 
ip{t^q) G C°^(]R X M^-^jC^ ). But then hnearity and the finite propagation speed 
(Proposition 1.1 in ^) imply that il){t,q) G C°°(M x RS^.C^'^) also for V^(0) G 
C°°(R3^,C^'^) n L2(R3^,C^'^). Hence, the Dirac current ^ satisfies (dH). Since 
I J| < cy/N j^, Corollary ^ implies the theorem. D 

Corollary 2 Let now Q = R'^ \ U^-^^Se, where Se is a hyperplane with codimension > 2 
for i = 1, . . . ,m, and let j = (j°, J) be a current as defined in (fTT|) . Suppose that there 
is a global bound c on velocities, \J\ < cj^, and that J and the first order derivatives of 
j° are bounded on bounded sets. 

Then for fiQ-alniost all q G R'^, the solution of (jH)) starting at Q{0) = q exists for all 
times, and the family of measures fit is equivariant. 

Proof. First note that R'^ \ Q is a Lebesgue-null set and hence also a /zo-null set. For 
g G Q we apply Theorem ^ The conditions (|15|) and (jl6|) of Theorem ^ follow as in the 
proof of Corollary ^ using the fact that J and the derivatives of j° are locally bounded. 
To check (fTTj) . let de be the dimension of Se and assume without loss of generality 
that Si contains the origin. Then with \J\ < C on B^^^j,, the ball of radius r + cT 
around the origin in R"^, and B'^^^j. c B^'^^j, x B^~Jp we find that 

^ dt [ dq y^!''^^}, <T f dx f dyfj <oo. 
J dist(g, S^) J J \y\ 






D 



5 Global Existence of Bohmian Mechanics 

We now apply Theorem^ to Bohmian mechanics and consider the abstract Hamiltonian 

Ho = -\ (m-^(V, - iA(g)))' 1^. + V{q) , D{Ho) = C^{Q,C^) , (22) 

where, for the moment, A G ifioc(R'^, R"') and V G Ll^^{Q,Y{.eim.{C^)) . The mass 
matrix m = diag(mi, . . . , m^) has positive entries mi > 0. These conditions assure 
that Hq is well defined and symmetric on C^{Q, C^). Since Hq commutes with complex 
conjugation, Hq has at least on self-adjoint extension. We also assume that Q = R"' \ 
U^-^Se where each Si is a (d — 3)-dimensional hyperplane in R'^. As to be explained 
in the example below, for d = 3N the coincidence set of A^ particles moving in R^ has 
exactly this structure and therefore singular pair-potentials like the Coulomb potential 
are included. In these abstract terms the Bohmian equation of motion reads 

§M.„-.,„,*:(^_^(,Q(,)). (23) 
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Theorem 3 Let H be a self-adjoint extension of Hq as in ()22|) with domain D{H). 
Suppose that for some ip{0) G D{H) with ||'?/'(0)|| = 1 the solution il){t) = e~'*^'?/'(0) of 
the Schrodinger equation satisfies 

(i) tpeC\Rx QX''), 

(ii) for every T > there is a constant Ct < C)0 such that 

I dt (lllVV^(t)lf +|||A^(t)|f +||A- WWf) <Ct. 

Then the solution Qq{t) of (j^ with Qq{Q) = q exists globally in time for almost all 
g e M"^ relative to the measure fio{dq) = \ip{0,q)\'^dq, and the \ipit)\'^ distributions are 
equivariant. 

Remark. 

9. Note that condition (i) in Theorem IHl is typically satisfied only if the potentials A 
and V are sufficiently smooth on Q, more than we required after (J22I)- We decided 
to state the condition in terms of ip since the exact type of smoothness required 
for A and V depends on, among other factors, the dimension d. 

Proof of Theorem First note that M"^ \ Q is a Lebesgue-null set and hence also a 
y^o-nuU set. For g G Q we apply Theorem [H According to Section H^ and by virtue of 
(i), the Schrodinger current 

jit,q) = {r{t,q)^{t,q), m-ilm^*(t, g)(V, - iA(g)) ^(t, g)) 

satisfies (fTT|l . We now check (fTH|) . (fTHI) and (fTTjl . in order to prove existence for positive 
times. For negative times one concludes analogously by applying exactly the same 
arguments to the time reversed current. 

With %l){t) = e~^^^tlj{0), the Cauchy-Schwarz inequality, and (ii) we obtain 

dq\dtj\t,q)\ = [ dq\dtr{t,qMt,q)\<2 [ dq\r{t,q)H^{t,q)\ 

Q JRd- jRd. 

< 2\\Hm\\ = nHi^m ■ 

For the second term in (|15|) we find, after a straightforward computation involving 
Cauchy-Schwarz first on C'^ and then on L'^{M.'^) and finally on L^([0,T]), that 



T T 

dt [dq l^.\/f{t,q) 



< — I dt I dq{\V^{t,q)\' + m,q)\\Aiq)-V^{t,q)\) 



Q\.Ai Rd 

Ct + \/TCt 
~ mo 

where m-o = min{mi, . . . , nid}- Hence, (11 5j] holds. Analogously p(i|l follows from 

T T 

dt fdq\J{t,q)\< [dt— [ dqm,q)\{\V^|J{t,q)\ + \A{q)^|J{t,q)\) <^^^ 
J J mo J mo 

Q Rd 
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We now come to (jTTj). Since S'^ is a (li — 3)-dimensional hyperplane, it can be written 
as Si = {q eR'^ : yi{q) = a/,} with ?/£ : M'^ ^ R^, q ^ {q ■ yj, q ■ yj, q ■ yf) where y\, yj, 
yf are 3 orthogonal unit vectors normal to the hyperplane S^ and a^ G M^ a constant. 
The distance to the hyperplane is given by dist(g, S() = \ye{q) — (ie\. 

To prove (fTTj) for 6 = oo, we use the generalized Hardy inequality introduced in |3], 
Equation (25). It states that for all (p ^ H^{M.'^, C), the first Sobolev space, 



Hence, 



[da . ^ff . <l i^WM,)?. 



Jq dist{q,Se) mo J J^d \ye{q) ~ ae\ 



< _Ljdt f ,jV^(t,g)|(|V^(t,g)| + |A(g)^(t,g)|) 



mo J 7]Rd \ye{q) 





T 



< — fdt(f dq y^'^^'' ) (|||V^(t)||| + |||A^(t)|||) 

mo J yJut \yi{q) -<^A 




T 



1 r nr^ 

< — \ dt{2\\ |V^(t)| f + II |VV^(t)| II II |A^(t)| II) < ^ 
mo J mo 





D 



We shall not try to verify the assumptions of Theorem IHl under as general as possible 
conditions on A and V. Instead we consider two examples where they can be checked 
without too much effort. 

Our first example concerns a molecular system in external fields. More precisely 
we consider A^ electrons in M.^ with configuration q = {q^, . . . , q^) G R^^ interacting 
through Coulomb potentials 

N~l N ^ 

in the electric potential 



i=i j=i+i i^« ^il 



Ku(g) = -EE 
i=i j=i 



N M 

^3 



\<l^-Zj\ 



of M static nuclei located at Zj G R? with charges Zj, j = 1, . . . , M. Furthermore we 
allow for an external magnetic field B{x) = V x A{x) with A G C°°(R^,M^) such that 
V ■ A = and B and A are bounded. The Hamiltonian of the system thus is 

(N \ ^ 
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with domain D(H^o\) = -ff^(ffi^^, (C^)*^^). Here cTi is the vector of Pauh matrices acting 
on the spin index of particle i. It is well known that V^i, V^^, and Vg are infinitesimally 
bounded with respect to A^. Hence -ffmoi = ~2^q + -^ with 

(N \ N 

-i 5^(2iA(qJ ■ V,^ - A{qf) + VM + Ku(g) l(c^)«iv - Y, B{q,) ■ cr, 

is self-adjoint by virtue of Kato's theorem. 

Corollary 3 Let ^(t) = e-"^-°i^(0) with ^(0) G C°°{H^oi) = r]^=iD{{HmoiT) and 
11^(0)11 = 1. Then the Bohmian trajectories Qq{t) exist globally in time for almost all 
q G M?^ relative to the measure |^(0, q)\'^dq, and the |^(t)p distributions are equivariant. 

Proof. First note that -ffmoi is of the form ()22|1 with d = 3N and /c = 2^. The configu- 
ration space of the system is 

Q = R'^\ ((uf^T^ Uf=„ {q G M^^: q, = q,}) U (uf=, ufU U ^ M^^: q. = z,})) , 

where the A^(A^ — l)/2 electron-electron and the NM electron-nucleus coincidence hy- 
perplanes are all (3A^ — 3)-dimensional. As remarked above, -ffmoi is self-adjoint on 
H'^[M.^^, (C^)*^^) and thus satisfies the hypotheses of Theorem El Hence it suffices to 
check that ip{t) satisfies the hypotheses (i) and (ii) of Theorem El As for (i), note that 
all potentials in (j24p are C°° on Q. Then methods of elliptic regularity can be applied 
to show that for ^(0) G C°°(ifmoi) the solution of the Schrodinger equation satisfies 
ip G C°°(R X Q). For details see the appendix in 4]. Finally notice that, since A is 
assumed to be bounded and since ||'?/'(i)|| = ||-?/'(0)||, (n) follows if we can show that the 
kinetic energy || |V^(t)| || remains bounded. This is also standard but we give the short 
argument anyway: since R is infinitesimally bounded with respect to A, there are con- 
stants < a < 1 and 6 > such that ||i?0|| < a|||A0|| +6||0|| for all 4> e H^ = D{H^oi)- 
Hence 

iiA^(t)ii = 2\\c^A-R+R)m\\<nH^m+nRm\\ 

< 2||i7^(t)||+a||AV^(t)||+26||V^(t)|| 
together with \\Hip{t)\\ = ||i/^(0)|| and ||V^(t)|| = ||^(0)|| implies 

IIAv,(i)||<^l|g^-(°)ll + ^''IN'(°)ll,c. 

1 — a 
But then also 

II |V^(t)| f = (V^(t),-V^(t)) = -(^(t), A^(t)) < ||^(t)|| ||A^(t)|| < 11^(0)11 C. 

D 

The last corollary coincides exactly with the result of |1] (see their Corollary 3.2). 
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Corollary 4 In ^ let k = 1, A = and V = Vi + V2 e C°°(Q,C), where Vi 
is hounded below and V2 "is —^A-form bounded with relative bound < 1. Then the 
form sum H = — |A + V is a self-adjoint extension of Hq and for ilj{t) = e~^*^'?/'(0) 
with V'(O) G C°^{H) = n^=iD(iJ"), 11^(0)11 = 1, the Bohmian trajectories Qg{t) exist 
globally in time for almost all q eM.'^ relative to the measure \ip{0, q)\'^dq, and the |^(t)P 
distributions are equivariant. 

Proof. For the statement about the form sum see [TSj. Again, as shown in the appendix of 
|3], elliptic regularity implies that ip G C°°(]R x Q). Hence, in order to apply Theorem El 
it suffices to show that that || |V^(t)| || remains bounded. This follows by an argument 
analogous to the one given in the proof of Corollary El For the details see the proof of 
Corollary 3.2 in m. D 
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